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1.  INTRODUCTION 


In  this  chapter,  some  partially  adaptive  space-time  processing  (PSTAP)  schemes  [1]  are 
compared  according  to  the  performance  obtained  using  a  measured  radar  data  set  called  multi¬ 
channel  airborne  measurements  (MCARM).  Algorithms  studied  include  adaptive  displaced  phase 
centered  antenna  (ADPCA),  factored  post  Doppler  (FTS),  extended  factored  approach  (EFA), 
joint-domain  localized  (JDL),  subarraying  ADPCA,  subarraying  EFA,  subarraying  FTS,  and 
beamspace  ADPCA.  We  use  normalized  test  statistics  which  achieve  constant  false  alarm 
probability  in  homogenous  clutter.  The  performance  of  the  algorithms  is  studied  in  section  2. 

1.1  Subarraying  ADPCA 


Let  N  be  the  number  of  antenna  elements,  M  be  the  number  of  pulses  per  coherent 
processing  interval  (CPI),  Ks  be  the  number  of  beams  formed  and,  Kt  be  the  number  of  pulses 
which  will  be  adaptively  processed. 

First  the  pre-processing  performs  the  beamforming  and  forms  groups  of  Kt  pulses  from 
Ks  beam  outputs  as 

Xk(p)=(Ap®G)HXk  ;  p  =  0,l,2,....,P-l  (1) 


where  Xk  is  the  space-time  snapshot  organized  as 


^k  1*1.1*  i  X2,1,k  » * "  »  XN,1,k  »  X1,2,kt  ' '  ■  XN,M,k  F 


(2) 


where  Xyjk  is  the  observation  from  the  i*  antenna  element,  j*  pluse  and  k*  range  cell  [2],  G  is 
the  NxKs  beamformer  matrix  ,  and 

0pxK, 

'k, 

°(M-Kt-p><Kt 


I  2 


(3) 


where  the  notation  0qxm  refers  to  an  qxm  matrix  of  zeros  and  IKt  is  a  KtxKt  identity  matrix.  It  is 
possible  to  form  sub-CPIs  which  overlap  by  any  number  of  pulses  as  in  [2].  In  this  report  we 
focus  on  the  case  where  Kt  is  set  to  3  and  consecutive  sub-CPIs  overlap  two  pulses.  We  also 
focus  on  Kg=3  in  this  report. 


The  beamformer  matrix  mentioned  in  [1]  (See  equation  (256)  on  page  14  in  [1].)  is  used  in 
our  processing.  This  matrix  is 
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where  g  =  (go,  gi,  g2, . . . ,  gN>-i)T  is  a  vector  of  scalar  components  and  N’=  N-Ks+1 .  The  vector 
g  could  be  defined  as  any  of  the  popular  windows  or  tapering  from  the  DSP  [3]  or  radar  [4] 
literature.  Here  we  take  g  as  either  a  uniform  window  (all  ones)  or  the  hamming  window  in  our 
tests.  The  uniform  window  gave  better  results  than  the  hamming  window  so  we  present  results  for 
the  uniform  window.  The  G  matrix  in  (4)  represents  a  case  where  the  beams  use  identical 
beamforming  using  a  length  N’  element  subaperture,  with  each  subaperture  shifted  by  one 
element  from  the  previous  one.  The  resulting  processing  is  sometimes  called  subarraying. 

Each  vector  produced  by  the  pre-processing  (1)  will  be  adaptively  processed  as 

yk(p)= sMR;’(p)xk(p)/<i)  (5) 


where 

S  =  (lKt®G)HS. 

(6) 

with  Seas 

s.  =S,®8. 

(7) 
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where  Ss  is  the  Nxl  spatial  target  vector  [1]  and  S,  is  a  K,xl  vector  of  binomial  coefficients  of 
altering  sign  normally  used  in  ADPCA  [2].  The  term  <£>  in  the  denominator  of  (5)  is  the 
normalization  to  provide  CFAR  in  homogeneous  clutter  and  is  given  by 

3>  =  S"Rk'S  (8) 

Rk  is  the  interference-plus-noise  covariance  matrix  estimated  from  the  pre-processed  data 
from  the  Q  adjacent  range  cells  (reference  data),  excluding  the  cell-  under-test  and  the  two  closest 
cells  on  each  side  of  the  cell-under-test,  as  follows 


a  k+Q/2+2  - 

rk(p)  =  ^  2  Xi(p)X,(p)H 

U  i=k-Q/2-2,i*k-2,k-1,k,k+1,k+2 


(9) 


After  the  adaptive  processing,  post  processing  is  performed  which  can  be  expressed  as 


Yk 


-iT 

yk(o)>yk(i) . yk(P-1) 


(10) 


= e  vk 


(id 


assuming  the  target  is  in  the  mth  Doppler  bin.  In  this  case  fm  is  the  m*11  column  of  a  PxP  Doppler 
filter  matrix  F,  and  Zk>m  is  the  final  output  whose  magnitude  will  be  compared  to  a  threshold  to 
produce  a  decision.  In  the  case  of  plain  ADPCA  (not  subarraying)  we  set  G  =  In  . 


1.2  Subarraying  FTS  (BFTS) 

In  subarraying  factored  post-Doppler  (FTS)  STAP,  the  pre-processing  is 
described  by  (1)  with  Ap=  fp  and  P=M.  fp  is  the  p‘h  column  of  an  MxM  Doppler  filter 
matrix  F=[f0,fi,f2> . 4m-i1  G  is  the  NxKs  beamforming  matrix  as  in  subarraying  ADPCA.  This 


4 


pre-processing  transforms  the  signal  into  Doppler  space.  To  test  for  a  target  at  a  particular 
Doppler  frequency,  only  one  pre-processing  from  (1)  needs  to  be  calculated.  In  FTS  a  single 
pulse  at  a  time  is  processed  (Kt=l),  and  p  indicates  the  Doppler  bin  in  question.  The  adaptive 
processing  in  (5)  is  employed  with  the  steering  vector  S  defined  as  in  (6)  with  Kt=l.  Here  Se  is 
set  to  Ss  in  (7),  since  Kt  is  1.  Post  processing  is  not  usually  employed,  so  |yk  (p)|  is  compared  to  a 
threshold  to  test  for  a  target  in  the  p,h  Doppler  bin.  In  the  case  of  plain  FTS  we  set  G  =  In  . 

1.3  Subarraying  EFA  (BEFA) 

The  subarraying  extended  factored  approach  (BEFA)  is  a  slight  extension  of  the 
subarraying  FTS  approach.  In  BEFA,  adaptive  processing  is  applied  to  several  adjacent  Doppler 
bins  instead  of  just  one.  In  our  subarraying  EFA  implementation,  the  scheme  adapts  over  3 
adjacent  bins.  The  pre-processing  is  as  in  (1)  with  Ap=[fp.i,  fp,  fp+1].  The  pre-processing 
performs  both  transformation  and  selection.  The  adaptive  processing  is  as  in  subarraying 
ADPCA  with  St=  [0  1  0] T.  Post  processing  is  not  employed.  In  the  case  of  plain  EFA  we  set  G  = 
In. 

1.4  Beamspace  ADPCA  (BeamAD) 

In  beamspace  ADPCA  ,  the  pre-processing  can  be  defined  as  in  (1)  with  Ap  as  in  (3),  and 
G=[  fN,i,  fN>2,  fNjKs  ],  where  fNJ  j=l,  K*  are  the  columns  of  NxN  DFT  matrix 
corresponding  to  Kg  angle  bins.  For  simplicity  assume  Ks=3  and  the  target  to  be  detected  is  in  the 
2nd  angle  bin.  The  adaptive  processing  is  as  in  subarraying  ADPCA  with  S  =  [0  1  0]®S,, 
where  St  is  a  Kt  x  1  vector  of  binomial  coefficients  of  altering  sign  as  in  subarraying  ADPCA.  In 
our  tests  Ks  and  Kt  were  set  to  3.  The  post  processing  is  as  in  subarraying  ADPCA  as  described 
by  (10)  and  (11). 

1.5  Joint-domaint  localized  approach  (JDL) 

In  JDL,  the  pre-processor  performs  two  dimensional  transformation  and  selection.  The 
data  is  transformed  from  the  space-time  domain  into  the  angle-Doppler  domain.  This  pre¬ 
processing  can  be  described  as  in  (1),  with  Ap=[  fm>1,  fm,2, ...»  fm,Kt],  where  fmJ  j=l, . .  K,  are 

the  columns  of  an  MxM  DFT  matrix  corresponding  to  the  Kt  Doppler  bins  and  with  G-[  fn.1,  f„4, 
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],  where  f„j  j=l,  Ks  are  the  columns  of  NxN  DFT  matrix  corresponding  to  the  Ks 
angle  bins.  As  for  FTS  and  EFA  only  the  post  processing  corresponding  to  a  single  p  must  be 
calculated  to  test  for  a  target  at  a  particular  Doppler  and  Angular  frequency. 

If  we  again  focus  on  Kt  =  Ks  =3  and  consider  the  case  where  the  target  to  be  detected  is 
in  2nd  angle  and  the  2nd  Doppler  bin  then,  the  adaptive  processing  can  be  described  as  in  (5)  with 

S  =  [0  1  0]®[0  1  0]  (12) 

More  precisely,  the  steering  vector  has  all  its  entries  equal  to  zero  except  for  the  one 
corresponding  to  the  angle  and  Doppler  bin  of  the  target.  No  post-processing  is  employed  for 
JDL. 


2.  REAL  DATA  PERFORMANCE 


To  test  the  STAP  algorithms  described  in  section  1,  we  use  data  that  comes  from  the 
MCARM  database  flight  5  acquisition  575.  See  [5]  and  [6]  for  detailed  information  about  the 
MCARM  program  and  the  data.  For  each  experiment,  a  single  target  signal  with  amplitude  0.05, 
a  particular  normalized  Doppler  frequency  and  a  particular  spatial  frequency  was  inserted  in  a 
particular  range  bin.  Reference  data  are  selected  from  consecutive  cells  on  each  side  of  the  cell- 
under-test  (see  Appendix  for  complete  details  on  all  schemes).  We  employ  normalized  test 
statistics,  which  provide  a  constant  false  alarm  rate  (CFAR)  characteristic  for  homogenous  clutter 
[7]. 


For  each  example,  we  provide  plots  of  the  magnitude  of  the  normalized  test  statistics  for 
a  set  of  ranges  including  the  target  range.  We  judge  a  scheme  by  how  large  the  test  statistic  is  at 
the  target  range  in  comparison  to  other  ranges. 

In  the  first  example,  we  inserted  a  target  at  range  bin  150  for  the  cases  shown  in  Table  1. 
The  location  of  the  targets  and  an  estimate  of  the  clutter  (plus  noise)  power  spectral  density  (psd) 
are  illustrated  in  Fig.  1.  As  visible  from  Fig.  1,  the  psd  estimate  used  is  rather  crude  and  is 
provided  to  give  a  rough  description  of  the  clutter  environment.  In  the  estimate  no  neighboring 
range  cells  are  averaged  and  a  Blackman  window  is  used.  Some  of  the  artifacts  can  be  removed 
by  averaging,  but  this  was  not  considered  necessary  in  this  case. 
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Case 

Normalized  Doppler 

Frequency 

Normalized  Spatial 

Frequency 

a 

-0.2656 

-0.2656 

b 

-0.0312 

-0.0312 

c 

0.125 

0.203 

d 

-0.0312 

-0.1875 

e 

-0.0312 

0.203 

f 

0.3593 

-0.1875 

g 

-0.1875 

0.3593 

Table  1 )  Test  cases  for  each  example. 

A  summary  of  the  results  is  given  in  Table  2.  Fig.  2  through  Fig.  8  presents  the  results  for 
most  of  the  schemes  tested.  FTS  and  subarraying  FTS  generally  perform  poorly,  so  their  results 
are  not  shown.  In  this  example,  generally  JDL  provides  best  results.  JDL  is  best  in  every  case 
except  cases  d  and  e.  In  case  d,  BeamAD  slightly  outperforms  JDL  and  ADPCA,  but  the 
difference  is  quite  small.  In  case  e,  EFA  outperforms  the  other  schemes,  however  JDL  also 
performs  well.  The  apparently  large  interference  indicated  by  Fig.  1  near  the  normalized  Doppler 
frequency  -0.0312  could  be  the  reason  why  JDL  is  not  best  in  these  two  cases.  These  are 
apparently  difficult  cases  where  no  scheme  can  really  excel.  Table  2  shows  that  case  b  is  also 
quite  difficult,  but  here  JDL  performs  much  better  than  the  other  schemes.  The  extra  clutter  ridges 
in  Fig.  1  are  discussed  in  [8]  and  [9]. 

In  the  second  example,  we  inserted  a  target  at  range  bin  350.  We  present  results  for  the 
same  cases  in  Table  L  The  best  three  schemes  for  all  the  cases  are  given  in  Table  3.  The 
normalized  test  statistics  for  the  six  best  schemes  tested  for  each  case  in  Table  1  are  given  in 
Figures  10  through  16.  The  results  indicate  that  none  of  the  schemes  always  outperforms  all  the 
others.  However,  post-Doppler  algorithms  are  generally  better  than  pre-Doppler  algorithms. 
Either  JDL  or  EFA  were  best  in  all  but  case  b.  In  case  b,  where  the  target  is  inserted  in  the  largest 
clutter  of  all  cases,  BEFA  is  only  slightly  better  than  EFA  and  JDL. 

Next,  we  inserted  a  target  at  range  bin  415.  We  present  results  for  the  same  cases  as  in 
previous  examples.  The  location  of  the  targets  and  an  estimate  of  the  clutter  psd  is  given  in  Fig. 
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17.  A  summary  of  the  results  is  provided  in  Table  4.  The  normalized  test  statistics  are  given  in 
Fig.  18  through  Fig.  24.  Here  for  cases  b,  c,  d,  f,  and  g  either  JDL,  or  BEFA  provide  best 
performance.  In  the  other  cases,  EFA  is  best  and  JDL  also  performs  well.  JDL  performs  well  in 
every  case  except  case  b.  Even  in  this  case,  its  performance  is  not  bad.  EFA  and  its  beamspace 
version  are  best  in  some  cases  and  near-best  in  others.  ADPCA  and  its  beamspace  version  give 
good  performance  in  number  of  cases,  but  these  schemes  were  never  best  in  this  example. 

Finally,  we  test  the  same  cases  as  before  when  the  target  is  inserted  in  range  bin  500. 
The  locations  of  the  targets  and  an  estimate  of  the  clutter  psd  is  given  in  Fig.  18.  A  summary  of 
the  results  is  given  in  Table  5.  Fig.  26  through  Fig.  32  provide  the  test  statistics  for  the  six  best 
schemes  in  each  case.  Again  the  beamspace  post-doppler  algorithms,  JDL  and  BEFA,  outperform 
the  element  space  algorithms  except  for  case  e  where  ADPCA  is  best. 


4.  CONCLUSIONS 

In  our  tests,  JDL  and  EFA  generally  perform  very  well.  Subarrayed  EFA  also  shows 
good  performance  in  many  situations.  The  common  element  these  schemes  share  is  post  doppler 
processing.  This  type  of  processing,  used  in  the  correct  way,  appears  to  be  superior  in  these  real 
data  cases.  This  appears  to  be  the  major  result  of  this  study.  To  a  much  lesser  extent,  schemes 
which  estimate  fewer  parameters  (lower  degree  of  freedom  (DOF)  schemes)  do  appear  to  provide 
some  benefits  in  non-homogeneous  cases.  The  number  of  parameters  estimated  for  each  scheme 
is  listed  in  Appendix  1.  Some  further  discussion  and  analysis  of  non-homogeneous  cases  is  given 
in  Appendix  2.  When  there  is  a  strong  interference  near  the  target,  as  in  case  d  and  e  in  Table  5 
for  example,  ADPCA  performs  well  and  sometimes  outperforms  all  the  other  schemes.  The 
reason  appears  to  be  related  to  the  extra  whitening  provided  by  its  steering  vector. 
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Figure  1)  Power  spectrum  plot  of  range  150 


Case 

The  3  best 

IEH 

Figure 

schemes 

M 

a 

JDL 

36 

-0.2656 

-0.2656 

BEFA 

25 

2 

BeamAD 

16 

■ 

-0.0312 

-0.0312 

JDL 

10 

3 

c 

JDL 

30 

0.125 

0.203 

EFA 

12 

4 

d 

BeamAD 

8 

-0.0312 

-0.1875 

JDL 

7 

5 

ADPCA 

6 

e 

EFA 

20 

-0.0312 

0.203 

JDL 

12 

6 

BeamAD 

5 

f 

JDL 

43 

0.3593 

-0.1875 

EFA 

38 

7 

BeamAD 

30 

9 

JDL 

38 

-0.1875 

0.3593 

BEFA 

18 

8 

BeamAD 

7 

Table  2)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range  150. 
(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range) 
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ADPCA 


BDPCA 


e)  EFA  f)BEFA 

Figure  2)  Performance  comparison  when  target  is  inserted  at  range150,  case  a. 


10 


e)  EFA  f)  BEFA 

Figure  3)  Performance  comparison  when  target  is  inserted  at  range  150,  case  b. 
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Normalized  spatial  frequecy 


Figure  9)  Power  spectrum  plot  of  range  350. 


Case 

The  3  best 
schemes 

D 

Figure 

EFA 

24 

a 

-0.2656 

-0.2656 

JDL 

18 

10 

■ 

BEFA 

5 

BEFA 

15 

-0.0312 

-0.0312 

EFA 

12 

11 

JDL 

10 

JDL 

34 

0.125 

0.203 

EFA 

25 

12 

BEFA 

13 

HHB 

EFA 

28 

-0.0312 

-0.1875 

JDL 

20 

13 

ADPCA 

12 

EFA 

22 

e 

-0.0312 

0.203 

JDL 

18 

14 

MmH 

ADPCA 

7 

mm 

EFA 

48 

gllB 

0.3593 

-0.1875 

JDL 

43 

15 

MBIBt 

BeamAD 

24 

JDL 

33 

9 

-0.1875 

0.3593 

BEFA 

18 

16 

BeamAD 

14 

Table  3)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range  350. 
(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range.) 
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e)  EFA  f)  BEFA 

Figure'll)  Performance  comparison  when  target  is  inserted  at  range  350,  case  b. 


19 


c)  JDL 


d)  BeamAD 


e)  EFA  f)  BEFA 

Figure  1 5)  Performance  comparison  when  target  is  inserted  at  range  350,  case  f. 


c)  JDL 


d)  BeamAD 


e)  EFA 


f)  BEFA 


Figure16)  Performance  comparison  when  target  is  inserted  at  range  350,  case  g 


Normalized  spatial  frequency 


Figure  1 7)  Power  spectrum  plot  of  range  41 5 


Case 

The  3  best 
schemes 

H 

Figure 

a 

EFA 

25 

-0.2656 

-0.2656 

JDL 

22 

18 

ADPCA 

12 

1 

BEFA 

25 

-0.0312 

-0.0312 

BeamAD 

18 

19 

EFA,  ADPCA 

15 

■ 

BEFA 

30 

8 

0.125 

0.203 

JDL 

23 

20 

BeamAD 

12 

d 

JDL,  BEFA 

20 

-0.0312 

-0.1875 

BeamAD 

18 

21 

EFA,  ADPCA 

16 

e 

EFA 

28 

-0.0312 

0.203 

JDL 

24 

22 

ADPCA 

15 

JDL 

45 

8 

0.3593 

-0.1875 

BeamAD 

35 

23 

BEFA 

32 

g 

JDL 

42 

-0.1875 

0.3593 

EFA 

20 

24 

BeamAD 

18 

Table  4)  The  three  best  schemes  for  ail  the  cases  when  the  target  is  inserted  at  range  415. 
(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range) 
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Figure*!  8)  Performance  comparison  when  target  is  inserted  at  range  415,  case  a, 


2' 


a)  ADPCA 


b)  BDPCA 


c)  JDL 


d)  BeamAD 


EFA 


BEFA 


e)  EFA  f)  BEFA 

Figure  20)  Performance  comparison  when  target  is  at  range  415,  case  c. 


a)  ADPCA 


b)  BDPCA 


c)  JDL 


d)  Beam  AD 


EFA 


BEFA 


e)  EFA 


f)  BEFA 


Figure  21)  Performance  comparison  when  target  is  inserted  at  target  range  415,  case  d 


c)JDL 


d)BeamAD 


e)EFA  f)BEFA 


Figure  24)  Performance  comparison  when  target  is  inserted  at  range  415,  case  g 
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Normalized  spatial  frequency 


Figure  25)  Power  spectrum  plot  of  range  500 


Case 

Normalized  Doppler 
Frequency 

The  3  best 
schemes 

D 

Figure 

a 

JDL 

32 

•  0.2656 

-0.2656 

EFA 

13 

26 

ADPCA 

7 

b 

JDL 

22 

-  0.0312 

-0.0312 

EFA 

20 

27 

BEFA 

15 

c 

JDL 

36 

0.125 

0.203 

BEFA 

12 

28 

7 

d 

JDL 

26 

-0.0312 

-0.1875 

EFA 

20 

29 

ADPCA 

10 

e 

ADPCA 

22 

-0.0312 

0.203 

JDL 

18 

30 

EFA 

15 

f 

JDL 

38 

0.3593 

-0.1875 

EFA,  BEFA 

28 

31 

BeamAD 

23 

g 

BEFA,  JDL 

33 

-0.1875 

0.3593 

BeamAD 

17 

32 

EFA 

15 

Table  5)  The  three  best  schemes  for  all  the  cases  when  the  target  is  inserted  at  range  500. 

(D  is  the  approximate  difference  between  the  normalized  test  statistic  at  the  target  peak  and  the 
largest  peak  in  the  normalized  test  statistic  at  some  other  range) 
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a)  ADPCA 


b)  BDPCA 


JOL 


bvamAD 


Figure  27)  Performance  comparison  when  target  is  inserted  at  range  500,  case  b 
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470  480 


900  510  820  830  540 


a)  ADPCA 


c)  JDL 


470  480  480  800  510 


b)  BDPCA 


70  480  480  500  510  520 


d)  BeamAD 


170  480  490  500  510  520  530  540 


d)  EFA 


f)BEFA 


Figure  28)  Performance  comparison  when  target  is  inserted  at  range  500,  case  c. 


Figure  29)  Performance  comparison  when  target  is  inserted  at  range  500,  case  d 


e)  EFA  f)  BEFA 

Figure  31)  Performance  comparison  when  target  is  inserted  at  range  500,  case  f. 
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APPENDIX  1 


Scheme 

N 

Kt 

Ks 

Ap 

G 

Q 

ADPCA 

128 

22 

3 

■ 

[  0px3,  I3,  0(l25-p)x3  ]T 

m 

66 

EFA 

128 

22 

3 

[fl28,Dp-1ifl28,Dpifl28,Dp+l] 

m 

66 

FTS 

128 

22 

1 

■ 

f|28,Dp 

m 

22 

Beamspace 

ADPCA 

128 

22 

3 

3 

[  0px3,  I3,  0(l25-p)x3  ]T 

See 

Note 

10 

Subarraying 

ADPCA 

128 

22 

3 

3 

[  0px3,  la,  0(i25-p)x3  ]T 

■ 

10 

Subarraying 

EFA 

128 

22 

3 

3 

[fl28,Dp-1  ,fl28,Dp.fl28,Dp+l] 

■ 

10 

Subarraying 

FTS 

128 

22 

3 

3 

f  1 28, Dp 

4 

JDL 

128 

22 

3 

3 

[fl28,Dp-1.fl28,Dp.fl28,Dp+l] 

See 

Note 

10 

Table  6)  Parameters  for  comparison  tests  in  this  report. 
Notes: 

Q  is  used  in  (9),  in  this  report. 


G  is  the  beamforming  matrix  which  is  formed  as  described  in  (13).  When  we  use  the 
uniform  g  vector,  the  G  (Ks=3  )  is  as  folows 
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1  0  0 
1  1  0 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
1  1  1 
0  1  1 
0  0  1 


Note:  MCARM  data  is  not  collected  by  a  uniformly  spaced  linear  antenna  array,  so  its 
beamforming  matrix  is  relatively  complicated.  This  G  matrix  used  in  beamspace  ADPCA 
and  JDL  can  be  obtained  by  using  measured  steering  vector  information  provided  with 
the  MCARM  database.  See  the  discussion  on  modSA  at  [5]. 
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APPENDIX  2 


For  real  data  it  is  difficult  to  evaluate  if  there  is  a  mismatch  between  the  statistics  for 
the  reference  data  and  the  statistics  for  the  data  from  the  cell-under-test.  However  by  considering 
energy  fluctuations  over  range,  it  is  possible  to  partially  evaluate  this.  When  the  mismatch  is  near 
the  target  location,  it  could  degrade  the  performance  dramatically.  To  evaluate  how  stationary  the 
reference  data  is  and  to  determine  how  the  statistics  of  the  reference  data  and  the  statistics  for  the 
data  from  the  cell-under-test  are  mismatched  near  the  target,  we  plot  energy  fluctuations  over  33 
range  cells  on  each  side  of  the  target.  For  each  range  we  plot  the  average  over  1 1  neighboring 
Doppler  bins  and  11  neighboring  angle  bins  near  the  target  Doppler  and  angle  bin.  We  use  5 
Doppler  and  angle  bins  on  each  side  of  the  target.  The  results  are  given  in  Fig.  33  through  Fig.  36. 

In  Fig.  33  through  Fig.  36,  in  all  plots  the  data  seems  to  be  non-stationary.  But  we 
judge  that  when  the  number  and  the  duration  of  the  humps  and  the  nulls  increases  in  these  plots, 
the  degree  of  non-stationary  of  the  data  also  increases.  Generally  when  the  data  is  highly  non- 
stationary  the  beamspace  post-doppler  schemes  JDL  and  BEFA  outperform  EFA  and  the  other 
element  space  schemes. 
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g) 


Figure  36)  Energy  plots  for  case  a,  b,  c,  d,  e,  f,  g  f  when  the  target  is  located  at  range  500. 
(Example  4) 
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